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^ ■ Abstract 

' A wireless packet network is considered in which each user transmits a stream of packets to its destination. 

^ i The transmit power of each user interferes with the transmission of all other users. A convex cost function of the 

Q completion times of the user packets are minimized by optimally allocating the users' transmission power subject 
to their respective power constraints. It is shown that, at all ranges of SINR, completion time minimization can 
OO ' be formulated as a convex optimization problem and hence can be efficiently solved. When channel knowledge is 

imperfect, robust power control is considered based on the channel fading distribution subject to outage probability 
constraints. The problem is shown to be convex when the fading distribution is log-concave in exponentiated channel 
power gains; e.g., when each user is under independent Rayleigh, Nakagami, or log-normal fading. 



I. Introduction 



Power control plays a crucial role in the operation of a wireless network, as it strives to provide 
> maximum benefits to the users within the confines of available resources. In particular, in a wireless 
network, each user's transmit power interferes with the transmission of all other users; therefore, power 
allocation has a significant impact on the quality-of-service (QoS) experienced by the network users. 
The user benefits derived from a given power allocation assignment can be characterized by different 
CN performance metrics. Traditionally, a common performance metric is a utility function of the rates of the 

So: 

Q , However, maximizing a concave utility function over the feasible rate region in a wireless network is not 
necessarily a convex optimization problem [l]-[5]. In the high signal to interference-plus-noise (SINR) 
regime, rate utility maximization can be approximately formulated as a convex optimization problem [1], 
and hence it can be efficiently solved. For voice or video applications that need to maintain at least a 
5^ moderate minimum rate, the high-SINR regime is often an appropriate assumption. On the other hand, in 
wireless sensor networks or delay-insensitive data applications, a user may wish to transmit at arbitrarily 
low rates (i.e., at low SINRs) depending on the channel conditions. At medium to low SINR, rate utility 
maximization are handled by iterative approximation methods [1], [2]. In this paper, we study a different 
network performance metric that is motivated by data applications. In particular, we consider the scenario 
where each user transmits a stream of packets to its destination, and we wish to minimize a convex 
cost function of the user packet completion times. We show that completion time minimization and the 
corresponding optimal power allocation can be formulated as a convex optimization problem at all ranges 
of SINR. 

Moreover, we consider imperfect channel knowledge due to channel fading, and formulate the mini- 
mization as a stochastic programming problem. The channel gains are modeled as random variables: an 
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(SPAWARSYSCEN), San Diego under Contract Nos. N66001-08-C-2013 and N66001-06-C-2020; and under ITMANET subcontract # 
18870740-37362-C issued by Stanford University and supported by the DARPA. Chris T. K. Ng is supported by a Croucher Foundation 
Fellowship. 



2 



outage event occurs when the SINR realization falls below the transmitter's SINR target. Robust power 
control is considered where each user is subject to an outage probability constraint. We show that for 
a wide class of commonly used channel fading distributions, e.g., Rayleigh, Nakagami, and log-normal, 
robust power control can be posed a convex optimization problem. 

Optimal power control in wireless networks is studied in [1]; it shows that maximizing concave rate 
utility functions can be formulated as geometric programming (GP) problems in the high-SINR regime, 
which are convex and hence their solutions can be computed efficiently. In the medium- to low-SINR 
regime, [6] describes an iterative approximation method to maximize a concave rate utility function 
by solving a series of GPs, and [2] proposes an algorithm to maximize a linear function of the rates by 
approximating the feasible SINR region by a series of polyblocks. In [3], sufficient conditions are presented 
for the convexity of the feasible QoS region, with optional constraints on the allocation of user power. 
Utility maximization through joint optimization of adaptive modulation, rate allocation, and power control 
is investigated in [5]. In [7] power control is studied in frequency-selective Gaussian interference channels. 
Outage probabilities corresponding to different fading distributions for network users and interferers are 
derived in [8]. For interference-limited wireless networks, optimal power control is considered in [9] under 
Rayleigh fading subject to outage probability constraints. 

The rest of the paper is organized as follows. Section |Il] describes the channel model and the optimization 
framework. In Section [nil minimizing a cost function of the completion times is posed as a convex 
optimization problem. Section |IV] explains how the completion time region and its corresponding rate 
region are related. Section IVl considers robust power control against imperfect channel knowledge subject 
to outage probability constraints. Conclusions are presented in Section IVTl 



II. System Model 

A. Wireless Channels 

Consider the scenario in which M users are communicating in a wireless packet network. Each user i 
consists of a transmitter i and a corresponding receiver i, where i = 1, . . . , M. Transmitter i wishes to 
send a stream of equal-length packets to receiver i, where each packet has Li bits. We assume the complex 
baseband additive white Gaussian noise (AWGN) channel model between the transmitters and receivers: 



= 2^ ,/H,,X, + Z„ (1) 
i=i 

where Fj is the observed signal at receiver i, Xj is the signal sent by transmitter j, Hij is the channel 
power gain from transmitter j to receiver i, and Zj is independent Gaussian noise with power Ni. We 
consider coherent detection: thus Hij may be assumed real as detection is unaffected by phase change 
in the received signal. In the subsequent sections, we consider different channel knowledge assumptions 
where Hi/s may represent known constants or random variables. Suppose transmitter i has transmit 
power constraint Pj. When transmitter i transmits at a power level of Pi < Pi, in the capacity limit, the 
transmission rate Ri (bits per channel use) achieved by user i is given by: 

P, = log(l + 5i), (2) 

where log is base 2, and Si is the signal to interference-plus-noise ratio (SINR) at receiver i. When the 
interference from the other transmitters is treated as noise, the SINR at receiver i is: 

We assume each user has a sufficient number of packets backlogged at the transmitter so that it is always 
transmitting, and we assume Li is sufficiently large to allow transmission close to the channel capacity 
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Fig. 1. An AWGN wireless packet network with AI = 2 users. 

limit. The completion time (number of channel uses) of the transmission of each of user i's packet is 
given by: 

Ti = Li/Ri. (4) 
For example, an AWGN wireless packet network with M = 2 users is illustrated in Fig. [U 

B. Completion Time Cost Function 

We consider the problem of minimizing a convex cost function of the completion times Ti, . . . ,Tjv/ 
by optimally choosing the users' transmission power subject to the power constraints: < Pj < Pj, 
i = 1, . . . , M. Let J(T) be the completion time cost function, where T = [Ti . . . Tm] , and other vectors 
are denoted similarly in this paper. We assume J(T) is jointly-convex in Ti, . . . ,Tm', the convexity 
penalizes overlong completion times. For example, the following completion time cost functions are 
convex [10]: 

J,(T) = T[i] + --- + Tm (5) 
Jp(T) = ((Tir + ---+(TMr)'^ p>l. (6) 

In ([5]), T[j] denotes the ith largest component of T. Thus the cost function Jr(T) is the sum of the r longest 
completion times. As special cases, r = 1 represents the maximum completion time: max{Ti, . . . ,Tm}, 
and r = M represents the sum of the completion times: Tj. In the cost function Jp(T) is the 
ip-norm of the user completion times. 

III. Packet Completion Time Minimization 

A. Perfect Channel Estimation 

We first consider the scenario where the channel power gains i/jj's can be accurately estimated and 
they are known by all users: 

= G,j, i = l,...,M, j = l,...,M, (7) 

where Gi/s are known constants. In this case, the minimization of the completion time cost function 
J(T) can be mathematically formulated as the following optimization problem: 

minimize J(T) (8) 
over Ti,Ri,Pi 
subject to 

> L,/Ri (9) 
P. < log! 1+ . ^ 1 (10) 



0<Pi<Pi, (11) 
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where i = 1, . . . , M, and the problem data Gij E IR+, Pi, Ni, E are given. The constraint (flOl) in 
the optimization problem is not convex [l]-[3], [5]. However, we show in the next section (l8])-(fTT]) can 
be transformed into a convex optimization problem, and hence its solution can be efficiently computed. 



B. Convex Optimization 

To formulate the the completion time minimization problem given in (l8l)- (fTTI) as a convex optimization 
problem, we first rewrite the constraints (fTOl) as: 

Next we apply the change of variables: 

^^ = ln^,, P^ = \nP„ t = l,...,M, (14) 

where In is the natural logarithm. The completion time minimization in (|8l)- (fTT)) then becomes: 

minimize J(T) (15) 
over Ti,Si, Pi 
subject to 

Li 

Ti > ^ (16) 

log(l + exp(S,)) 



In 



|exp(^, -Pi-\n Gii + In N{) + ^ exp(5i -Pi + P, - In Gu + In G,,) | < (17) 

P^<\nPi, (18) 
where i = 1, . . . , M. Note that the constraint in (fTTl) follows from rewriting (fTSl) as: 

S.Pr^GT^Ni + S^Pr^PfiT^G,, < 1, (19) 

and taking logarithm on both sides after applying (fT4l) . The change of variables is similar to the transfor- 
mation techniques in geometric programming (GP) problems [1]. In particular, the log-sum-exp function 
in constraint (flTl) is convex [10]. The convexity of (fT6l) can be verified from its second-order conditions, 
as the right-hand side of (fT6l) is twice-differentiable and its second derivative is positive: 



(f r , ^,.-1 eMn2 2e^-ln(l + e^)i 

— log 1 + en = ^ ^ (20) 

dx^^ ^ (l + e-)2[ln(l + e-)]' 

> 0, (21) 

which follows from the inequality y > ln(l + y) for y > 0. Note that the transformation in (fT4l) does 
impose a slight loss of generality as we assume Pi ^ 0. Nevertheless, the formulation in (fT5l)-(fT8l) is 
otherwise valid for all ranges of SINR, and its solution can be efficiently computed by standard numerical 
methods in convex optimization [10]. Note that we may consider additional linear or convex constraints 
on T, and sum power constraints on subsets of P: they can be readily incorporated in the optimization 
problem without violating its convexity. 
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C. Fading Channels and Power Adaptation 

In this section we consider fading channels, i.e., the channel gains i/j/s in ([T]) experience random 
variations. In particular, we assume the channel gains can be characterized by a set of s G {1, ... , S} 
discrete fading states: 

{G*^^^ with probability pi g 
\ , Ep^ = 1' (22) 

G^-^) with probability ps ''"^ 

where H G M^^^^'^ = [Hij\ is the channel power gain matrix, and G G M|^^^^ = [G^j\ is defined 
similarly. For example, the discrete states may represent a finite set of quantized channel estimates. We 
consider slow fading where the duration of a fading state is long compared to the packet completion 
times. We assume the channel state s can be accurately estimated and it is known by all users. 

We first consider the case where each user can adapt its transmission power level according to the 
fading state. Suppose user i transmits at power level in fading state s, subject to the average power 
constraints: 

s 

E[P^]=$^P.P/'^<P^, z = l,...,M. (23) 

s=l 

To minimize a cost function of the expected completion times, the optimization problem can be formulated 



as: 



minimize J(E[ri], . . . , E[Tm]) (24) 
over E[T,],T(^\^f\p(^) 
subject to 

s 

E[T,] = ^P,t/^) (25) 

s=l 

tI'^ > , , (26) 

"log(l + ^f^) 

f^{s) p{s) 

< P/'^ (28) 
s 

J2PsPt^<P^, (29) 
s=l 

where s = 1, . . . , S; i = 1, . . . , M. The optimization in (l24l)-(|29l) can then be transformed into a convex 
optimization problem by similar techniques as described in Section IIII-BI Note that to minimize the 
expected value of the cost function, it can be handled similarly by replacing the objective function in (|24)) 
by: 

5 

E[J{T)] = J2PsJ{Tt\...,Tl^), (30) 
where convexity is preserved in the nonnegative weighted sum of convex functions. 
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In the case where the user cannot adapt its transmission power level to the fading state (i.e., the 
transmitter is under a short-term power constraint), the optimization problem is similar to (|24l) -(|29l). but 
with the average power constraint in (|29l) replaced by separate power constraints for each fading state: 

P^<hi s = l,...,S, t = l,...,M. (31) 

Note that under the short-term power constraints of (|3T1) . minimizing the expected cost function (|30l ) 
decomposes into S independent optimization problems: i.e., each of J{T^'^\ . . . ,Tj^f^), for s = 1, . . . ,S, 
can be minimized separately. 



IV. Relations to Rate Region 

In general, in a wireless network as defined in ©-(SI), minimizing a convex cost function J(T) of the 
completion times Ti, . . . , Tm is not equivalent to maximizing a concave utility function f/(R) of the rates 
Ri, . . . , Rm- In particular, maximizing f/(R) over the rate region is in general non-convex [3]: at high 
SINR it can be approximately formulated as a GP, and in the medium- to low-SINR regime there are 
iterative approximation methods [1], [2]. Suppose the cost function J+(T) is convex and nondecreasing in 
each argument Ti, then completion time minimization is a special case of rate utility maximization where 
the optimization problem can be formulated as convex. To see that minimizing J+iT) can be posed as a 
rate utility maximization problem, we define the corresponding rate utility function: 

UtCR) ^ -MLjRu . . . , Lm/Rm). (32) 

Note that minimizing J+{T) is equivalent to maximizing Ut(R), and the utility function Ut(R) is concave 
in R as prescribed by the convexity composition rules [10]. 

Nevertheless, in the converse, some rate utility maximization problems can be formulated as minimizing 
convex functions of the completion times. Suppose we minimize a positively weighted sum of the 
completion times: 

J^(T) =wiTi + --- + wmTm, who, l^w = 1, (33) 
where >z denotes component-wise inequality, then it is equivalent to maximizing: 



w', w'j. 



U,{n) = ^, (34) 

where Ud(R) is the utility function that corresponds to minimum potential delay fairness [11], with 
w'- = WiLi, i = 1, . . . , M. In addition, minimizing JwiT) in (|33]) is also equivalent to maximizing the 
weighted harmonic mean of the rates: 

t/.(R)=(|^ + --- + |^)" (35) 

Note that by applying Jensen's inequality on the convex function 1/x for x E we have 

I < !fl + . . . + (36) 

w[Ri + ■ ■ ■ + w'j^fRM Ri Rm 

which implies: 

w[Ri + --- + w'j,jRM>Uh{n). (37) 

Hence maximizing Uh(R) provides a lower bound to max w[Ri + ■ ■ ■ + w'j^,jRm, which represents a 
weighted throughput of the wireless network. In particular, the bound is tight when Ri = ■ ■ ■ = Rm as 
equality is achieved in (|36l) : therefore, maximizing the minimum rate in (|38]) below can be formulated as 
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Fig. 2. Completion time region (Pi = P2 = dB). 



a minimization of the convex cost function Jx{T) as given in (l39l) . which corresponds to the maximum 
completion time: 



f/„(R) = min{i?i, . . . , Rm} 
J^(T) = max{Ti, . . . , Tm}. 



(38) 
(39) 



Moreover, the entire rate region achievable under (fTOl)-(fTTI) can be characterized in terms of the 
corresponding completion time region. Specifically, the completion time region as characterized in (fT6l)- 
(fTSl) is convex, and its boundary are given by the minimizer of JwiT) in (l33l) over all vi^ ^ 0, l^w = 1. In 
turn, from the monotonicity of dH), each minimal completion time vector T* = argmin Jw{T^) corresponds 
to a maximal rate vector R* on the boundary on the rate region (fTOl)- (fTT)) . with R* = Li/T*, for 
i = 1, . . . , M. As a numerical example, we consider the following 2-user AWGN wireless packet network: 



G 



0.42 
0.63 



0.89 
0.15 



dB, 



N 



10 
10 

0' 




dB, 



(40) 



(41) 



with maximum completion time constraints: Tj < lOOLj. The completion time region and its corresponding 
rate region are shown in Fig. [2] and Fig. |3l respectively. Note that the power constraints as given in (|4TI) 
belong to the low SINR regime, where the high-SINR GP approximation does not readily apply. The 
completion time region is convex in Fig. |2l however, note that its rate region counterpart is non-convex 
as can be observed in Fig. |3l 



V. Robust Power Control 

A. Outage Probability Constraints 

In Section Uni we assume that the channel power gains Ht/s can be accurately estimated. However, 
in a fading environment where the channel estimates are updated not as fast as the channels vary, the 
transmitters may not know the -ffj/s perfectly. In this section, we consider the scenario where the channel 
gains are modeled as random variables: 

H = W ^ [W,,] , (42) 

where the users know the joint probability distribution of Wi/s but do not know their realization. As in 
Section ITlI-C[ we assume the duration of a fading state is long compared to the packet completion times. 
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Fig. 3. Rate region (Pi = P2 = OdB). 



Based on the channel distribution, each user i chooses a target SINR Si. Should the realized channel 
SINR fall below the target Si, the receiver cannot decode the transmitter's message, and it results in 
an outage event. To ensure the network operates with an acceptable level of reliability, we extend the 
completion time minimization framework in Section [ni] to additionally consider constraints imposed on 
the permissible probability of outage. Specifically, we minimize the completion time cost function subject 
to a set of outage probability constraints: G [0, 1], z = 1, . . . , M, where we stipulate that the probability 
of user i's transmission in outage not exceed g^. 

Incorporating the outage probability constraints gj's, the minimization of the completion time cost 
function is described by the following stochastic programming [12] problem: 

minimize J(T) (43) 
over Ti.Si, Pi 
subject to 

T. > (44) 
log(l + Si) 

Pr| .. , <sA<q. (45) 



0<Pi< Pi, (46) 



where i = 1, . . . , M. We show that the minimization in (I43l)-(I461) can be posed as a convex optimization 
problem for a wide class of channel fading distributions commonly considered in wireless communications. 

B. Robust Convex Optimization 

In terms of the transformed variables in (fT4l) . we first define the reliability function as: 

^i{Si, P) = Pr{user i not in outage} (47) 



PiS^W,, > ln|exp(5i -Pi + InNi) + ^exp(S, - + P, + 



(48) 



where (|48l) follows from rearranging (|45l) with the additional change of variables: 

%^lnW,„ t = l,...,M, j = l,...,M. (49) 
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Next, we characterize the reliability function in terms of the channel distribution. Let Wj be an M- 
component non-negative random vector that corresponds to the ith row of the channel matrix W: 

W,^ [Wa ... W,Mf, (50) 

and Wj G M.^ a realization of Wj. Under (|49l ). the transformed vectors Wj, w, are defined similarly. Let 
/wi(wj) denote the joint probability distribution function (pdf) of Wj. Theorem [T] below describes the 
sufficient condition that establishs the log-concavity of $i(S'j, P), under which (I43l)-(l46l) can be posed as 
the following convex optimization problem: 

minimize J(T) (51) 

over Ti,Si,Pi 

subject to 

Li 

Ti > z ^ (52) 

log(l + exp(5j)) 

ln$,(S„P)>ln(l-g,) (53) 

Pi<\nP,, (54) 

where i = 1, . . . , M. Let exp(w,;) denote component-wise exponentiation of the vector Wj. 

Theorem 1: The reliability function $j(S'i,P) is log-concave in Si,P if /w, (exp(wj)) is log-concave 
in Wi. 

The proof is given in the appendix. For example, the above condition is satisfied when each user expe- 
riences independent fading distributed as [13]: i) Rayleigh (richly scattered environments), ii) Nakagami 
(significant line-of-sight propagation), or iii) log-normal (shadowing due to signal attenuation). Therefore, 
in all these cases, completion time minimization subject to outage probability constraints can be formulated 
as convex optimization problem (|5T|)-(l54l). 

VL Conclusions 

In this paper we consider minimizing a convex function of the completion times of user packets by 
optimally allocating transmission power in a wireless network. We first focus on the scenario where the 
channel gains can be estimated accurately and are known by all users. We show that completion time 
minimization can be formulated as a convex optimization problem, and hence the corresponding optimal 
power allocation can be efficiently computed. The optimization formulation is valid for all ranges of SINR, 
which is especially pertinent for wireless sensor networks or delay-insensitive data applications where the 
users may transmit at moderate or low SINRs. Under fading channels with transmission power adaptation 
across fading states, an average power constraint can be incorporated into the optimization problem. We 
show that completion time minimization is a special case of rate utility maximization for which the 
optimization problem can be posed as convex. In particular, in a wireless network, although the feasible 
rate region is non-convex, the corresponding completion time region is shown to be convex. Finally, we 
consider robust power control under imperfect channel knowledge in fading channels. Completion times 
are minimized subject to outage probability constraints over the fading distribution, and we show that for 
a wide class of commonly used fading distributions, e.g., Rayleigh, Nakagami, and log-normal, robust 
power control can be posed as a convex optimization problem. 

Appendix 

Proof of Theorem [7] 

We first introduce, in terms of the transformed variables in (|49l) . the notation of W_j G R*^~^ 
representing the interfering channel random vector: 

W_, ^ [Wa . . . . . . WaiY , (55) 
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and w_i is a realization of W_j. The pdf of W„j is given by the marginal: 

/w ,(w^0 = / /w,(w.) dWu. (56) 
Next, conditioning on W_j = w_j, the reliability function in (l48l) is given by: 

MS^, P)= [ US^., P, w_,) dw_,. (57) 

For notational convenience, we write P, w_i) as a composition of: 




where F-^y^Xwu) in (|59l ) is the cumulative distribution function (cdf) of Wjj, and (wjj) is referred to 
as its complementary cdf, which is a nonincreasing function in Wa. 

Proof: In the construction of $i(S'i, P) in (l57l) . with the application of Lemma [T] below, log-concavity 
is preserved [10], [12] in the integration in (l57l) . (l56l) : multiplication in (1581) : complementary cdf in (|59| ): 
and composition of a logarithmically concave, nonincreasing function with a convex function in (|60l) . ■ 
Lemma 1: Under transformation (|49l ). /wi(wt) is log-concave in Wj. 

Proof: Consider the logarithm of the pdf of transformed random vector Wj: 

ln/w,(wi) = ln< exp(^^Wijj/w,(exp(wi)) I (61) 

j=i } 

M 

= ^Wij + \Yi fw, (exp(wi)) , (62) 
where log-concavity of /w, (exp(wj)) in Wj follows from assumption from Theorem [T] ■ 
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